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Motivated by recent experiments, we present a theoretical investigation of how the electro-osmotic 
flow occurring in a capillary is modified when its charged surfaces are coated by charged polymers. 
The theoretical treatment is based on a three dimensional model consisting of a ternary fluid-mixture, 
representing the solvent and two species for the ions, confined between two parallel charged plates 
decorated by a fixed array of scatterers representing the polymer coating. The electro-osmotic 
flow, generated by a constant electric field applied in a direction parallel to the plates, is studied 
numerically by means of Lattice Boltzmann simulations. In order to gain further understanding 
we performed a simple theoretical analysis by extending the Stokes-Smoluchowski equation to take 
into account the porosity induced by the polymers in the region adjacent the walls. We discuss 
the nature of the velocity profiles by focusing on the competing effects of the polymer charges and 
the frictional forces they exert. We show evidence of the flow reduction and of the flow inversion 
phenomenon when the polymer charge is opposite to the surface charge. By using the density of 
polymers and the surface charge as control variables, we propose a phase diagram that discriminates 
the direct and the reversed flow regimes and determine its dependence on the ionic concentration. 


I. INTRODUCTION 

Electrokinetic phenomena of fluids under conditions of extreme confinement are important to micro and nanofluidics 
and have a variety of applications ranging from fabrication of efficient nanotechnological devices, electrochemical 
energy storage, electrokinetic energy conversion, up to biomedical applications, such as separation and analysis of 
biological molecules and molecule delivery and sensing [IHS]. What makes micro-sized channels attractive is the 
large surface to volume ratio, so that the surface has a greater impact and some new phenomena arise, opening 
the possibility to develop new fluidic functionalities, since decreasing the scales increases the sensitivity of analytic 
techniques which are used in Lab on a chip (LOG) devices [ll[8]. 

The present paper investigates the effect of coating the inner walls by polymers on the electro-osmotic flows (EOF). 
In standard electroosmosis the motion of ions and their surrounding water molecules, induced by an applied electric 
field parallel to the charged walls of a capillary, generates a flow that extends outside the Debye electric double 
layer (EDL). The presence of a non uniform velocity profile associated with the EOF may represent a problem 
in capillary electrophoresis or in microfluidic devices used in protein analysis, because it increases dispersion and 
reduces resolution. On the contrary, in capillary electrochromatography the flow must be enhanced to produce high 
throughputs. It has been demonstrated by experiments, computer simulations pun] and theoretical arguments m 
that a modification of the chemical composition of the EDL by coating the walls with polymers leads to a consistent 
reduction of the mass flow and even to a reversal of the electroosmotic current P . This modification of the flow may 
occur for two reasons: the first is the drag force exerted by the polymer beads on the electrolyte solution and the 
second is the electric field originated by the polymer charges. Danger and coworkers [T2| observed that decreasing 
the EOF optimizes the resolution and analysis time in the separation of peptide mixtures in coated capillaries and 
investigated which polymers were more efficient to realize such a situation. 

In spite of the large amount of experimental work on the mechanisms controlling the EOF in polymer coated 

capillaries, a complete understanding of the problem is still missing m- The scaling theory of Harden et al. [18] 
focused on the interplay between the deformation of the adsorbed polyelectrolytes and the fluid motion. A series 
of Molecular dynamics (MD) simulations gave support to their predictions regarding the dependence of the polymer 
length on the coatings. The influence of the polymeric structure and of the solvent on the EOF was also investigated 
at length scales smaller than the Debye length [I9H^ . Slater and coworkers [22| [23] found by MD that the direction 
of the EOF can reverse in the region near the walls. 

We propose here an alternative simulation approach to determine the EOF in coated capillaries, based on the Lattice 
Boltzmann method (LBM) [21] and is supported by a straightforward theoretical treatment of the current modulation. 
We believe that our numerical approach can be generally employed for this type of systems and is complementary to 
particle-based simulations, the latter being computationally more expensive such that only a limited range of length 
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scales, electrolyte concentrations and polymer coatings can be explored. In fact, MD provides information at the 
molecular level at the price of the heavy computational effort required to track individual molecules and, in order to 
reach statistical accuracy for electrolytic solutions where the ionic densities are several orders of magnitude smaller 
than solvent density, demands for expensive ensemble averaging procedures |25j . 

Our modeling of the polymer coating is idealized and neglects the fact that polymers with one end fixed to the 
surface tend to deform due to the motion of the fluid. The polymer layer is assimilated to a region populated by 
an assembly of fixed random charged obstacles exerting a drag force on the fluid. A standard description of such a 
situation is represented by the phenomenological Brinkman equation [26) . often used to model the behavior of the 
fluid velocity, u, through highly permeable porous media as: 

? 7 e//V^u = Vp+|u (1) 

where Vp is the pressure gradient, rj and rjeff are the standard and effective dynamic viscosity, respectively, k is the 
permeabilty of the medium, with ? 7 e// = 0 corresponding to the Darcy equation. 

In the following we will determine the properties of the EOF in slit-like capillaries with coating represented at 
finer level with respect to the description of eq. by using the LBM approach: a formula equivalent to eq. will 
be obtained from the one-particle phase space distributions of the individual species as basic variables through an 
averaging procedure. The fluid is described as a a ternary mixture comprising a solvent species and two ionic species 
of opposite charge, as recently done in ref. m- The LBM provides a statistical description of the system under study, 
allows to compute density profiles and flow fields and embodies correctly the hydrodynamics and the electrostatics. 
To this purpose, we consider the simplest version of the LBM for electrolytes which is an extension of the Bhatnagar- 
Gross-Krook (BGK) relaxation-time model. It neglects the short range structure of the fluid and the non ideal 
gas features of the fluid besides the long range electrostatic interactions, but also displays undeniable advantages such 
as the simplicity of programming, speed of calculation and allowing to explore arbitrary channel geometries. 

This paper is organised as follows: in [IT] we illustrate the Boltzmann-like transport approach describing a ternary 
mixture comprising a solvent species and two ionic species of opposite charge. The non electrostatic part of the 
collision terms is treated within a single relaxation time approximation of the BGK type [5H]. We derive the balance 
equations for the momenta and the densities and recover the structure of equations for the charge, the mass and the 
momentum predicted by the macroscopic models. In |III| in order to gain further insight we perform a reduction of the 
Boltzmann coupled equations to a Stokes-Smoluchowski formulation and obtain the electric potential profile within 
the pore in the Debye-Huckel approximation [5] and show analytically how the mass current can be modulated by 
changing the properties of the polymer adsorbed on the surfaces. We study the resulting velocity field as a function 
of the relevant parameters. In sections [TV] and [V| we validate such results by solving numerically the LBM equations 
and treat the electric potential using the Poisson equation. We compare the numerical and theoretical results for 
the velocity profile and the potential and show that in the limit of small surface charge densities the agreement 
regarding the global properties of the system are satisfactory, whereas for larger densities the deviation is due to the 
non-linear nature of the self-consistent electric potential. Finally, in section [VT| we make some conclusive remarks and 
considerations. 


II. THEORETICAL BACKGROUND 


Let us consider a ternary mixture and denote the components by the superscript a — 0,where 0 identifies 
the solvent, — and -I- indicate the anions and the cations, respectively. The three types of particles have masses m“, 
charges z°‘e (with = 0), expressed in units of the electronic charge e, and are free to move in a slit-like channel 
whose parallel walls coincide with the planes at z = ±w, are impenetrable and carry a fixed surface charge of density 
per unit area E, which is negative in the cases here considered. 

Our goal is to determine the relevant dynamical observables of the system, such as the partial number density 
n“(r,t) and local velocity u“(r,t) starting from the one-particle phase space distributions /“(r,v,t) by performing 
the following projections in velocity space: 


( n“(r,t) 

^ n“(r,<)u“(r, t) 



/“(r,v,t). 


( 2 ) 


From these quantities one derives the local charge distribution: 


= e^z“n“(r,t), 

a 


( 3 ) 
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the global number density: 




( 4 ) 


and the average barycentric velocity of the fluid: 

u(r,t) = 




( 5 ) 


Hereafter, in the notation we will drop the dependence (r, t) for these fields for the sake of clarity and also assume 
the temperature to be constant. 

We now consider the evolution equation of the /“’s (see reference [55] for details) : 

0 d v) 

-r(r,v,t) + V.vr(r,v,t) + —(^^.r(r,v,t)) 
at av m“ 

= -w[/“(r,v,f) (r,v,t)] 

+ ^V<^(r).|-r(r,v,f). 

m“ av 

( 6 ) 

The left hand side of eq. [^represents the streaming contribution to the evolution, while the right hand side contains 
a BGK relaxation term representing the effect of the molecular collisions tending to restore on time scale the 
local equilibrium distribution represented by: 


?/’“(r,v,t) = n“[ 


2TTkBT 


exp^- 


^(v — u)" 
2kBT 


( 7 ) 


with ks is the Boltzmann constant. 

The last term in the r.h.s. of eq. describes the electrostatic coupling among the ions, treated within the Vlasov 
mean-field approximation, where one assumes that the electrostatic potential (j) is the solution of the Poisson equation: 




( 8 ) 


where ptot is the sum of the mobile ionic charge density p and the density due to fixed charges p/(r) distributed 
within the system. The confining walls impose the vanishing of the densities n“ for \z\ > w and of the velocities u“ 
at z = itw according to the no-slip boundary condition. The walls also contribute to the electric potential cj), due to 
their surface charges. The fixed surface charges appear through the Neumann boundary conditions on the gradient 
of the electrostatic potential: 




reS(r) 


S(r) 


( 9 ) 


where n is the local normal to the surface, 5'(r). The associated value of the potential at the wall is: 

^lreS(r) — C- 


( 10 ) 


We model the interaction between the polymers and the fluid by taking into account that the polymers are adsorbed 
irreversibly on the surface and exert two kinds of forces on the fluid: a) a drag force in the region where the mobile 
particles hit the polymer and their momenta are reduced; b) an electric force, because the polymer charges interact 
with the ions of the solution and modify the structure of the EDL. The polymer coating is idealized as a region of 
thickness 6 = {w — h) and volume Vs adjacent to the wall and containing Ns fixed charges described by the density 
held 


Ns 

E- 


P/(r) = ez, V(5(rfc -r) 


( 11 ) 


where r*, are the positions of the charges and Zs their valences, and act as sources of the potential via eq. The 
drag force is assumed to be proportional to the velocity of the moving particles and to the density of obstacles: 


drag 


(gv) 


= —k 


Ns 

pE 


v5{rk - r) 


( 12 ) 
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where kp is a measure of its strength, eq. [^reproduces the correct hydrodynamic behavior, in fact, by integrating 
w.r.t. the velocity in eq. [^we obtain the conservation law for the particle number of each species: 




+ V 


(n“u)+V • (n“(u“ - u))= 0, 


(13) 


where the last term in eq. is the so-called dissipative diffusion current, measuring the drift of the a-component 
with respect to the barycentric velocity. After multiplying by v and integrating eq. [^w.r.t. v, we obtain the balance 
equation for the density of momentum of the species a: 


- n“« - - m,)) 

Ns 

= ^ S{r - rk)u^n^ 

k^l 

—ez°‘n°‘V j(f> — — uj) 


(14) 


where 


y “ Uz)ivj - V, t) 


(15) 


represents the kinetic contribution of component a to the pressure tensor m- In eq. [15 and in the following, the 
convention of summing over repeated indices is assumed. 

The balance of the total moment of the fluid follows immediately by summing eq. |14|over the three species: 


Pm^i^i^j ^i^ij 
Ns 

PmkpUj ^ (5(r - rfc) - ^ ez^n^Vjcj) = 0 

k—1 ± 

(16) 


where Pm = is the local mass density and Vij = X)a i® total (kinetic) pressure. By a standard 

analysis m for species of equal masses {mp = m), one obtains the mutual diffusion coefficient D = the dynamic 
shear viscosity rj = and the total kinetic pressure: 


/ dui duj 


N-kBTn5,j v[g^+g 


2 duk^ \ 

Sdxk^J 


(17) 


III. ANALYTIC TREATMENT 

The presence of charges on the walls and on the fixed polymers gives rise to the formation of an EDL along the 
z direction. Under the action of an electric field parallel to the walls induced by the presence of two electrodes, the 
excess positive charge in the EDL generates a flow of positive ions towards the cathode and their motion is transmitted 
by shear forces to the rest of the fluid in the capillary. 

Before embarking on the full numerical approach to the problem, it is useful to study the solution of the electro¬ 
hydrodynamic equations eq. and eq. by using the linearized Debye-Huckel theory to treat the EDL and the 
Stokes approximation to describe the creeping flow regime, commonly realized in micro and nanofluidics set-ups. 
In this limit, the non-linear terms in the velocity in eq. can be neglected being much smaller than the viscous 
term and will be dropped in the rest of the paper. To proceed analytically we average over the distribution of 
charged obstacles, assuming that they are uniform within the two layers each of volume 14 where they are located: 
^ — r))- Such an averaging leads to smooth densities and velocity profiles compared to the noisy 

profiles observed in the simulations, but allows to exploit the statistical invariance with respect to translations in 
directions parallel to the walls. 
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The resulting stationary state of the total momentum eq. |16[ averaged over the polymer distribution, is a solution 
of the following equation for the fluid velocity: 




^xP - ez^n^ + jpm'diz)ux = 0 

± 


(18) 


where we introduced the friction constant 7 = and the driving electric field = —'S/x'<p{x) parallel to the 

walls. The function Q{z) = {9{w — h — z) + 9{z — w + h)), constructed as the sum of two Heaviside distributions, 
describes the presence of the polymers in two slabs adjacent the walls and WxP indicates the term ksTVxn associated 
with a possible pressure gradient along the x direction. We also separate the pressure tensor into its dissipative and 
non dissipative contributions by using eq. |17[ Under steady conditions the density distributions are constant along 
planes parallel to the walls and, neglecting steric effects due to the short range repulsive forces between molecules, 
can be determined with the help of eq. and by assuming the local equilibrium condition. One obtains the following 
self-consistent Poisson-Boltzmann equation for the electrostatic potential (j) between the plates: 


d?(j) 


eub . , 
= -smh 


ksT 


eZsUs 


e. 


i2 1 SZsTlg 

kjjCj) 0 , 


where nt is the density of the single ionic species in the bulk and k^ is the inverse Debye length: 


(19) 


Xd — kjj 



( 20 ) 


Notice that in the r.h.s. of eq. [^we separated the source into two contributions, the first taking into account the 
mobile charges and the second the charge of the obstacles. In the analytic treatment we employ the linear Debye- 
Huckel approximation for the electrostatic potential corresponding to the second approximated equality in eq. [T^ 
but in the numerical LBM treatment we consider the fully non-linear Poisson-Boltzmann equation. The solution of 
eq. with the properties of being continuous together with its first derivative at z = zLh and corresponding to the 
surface charge density E at z = iru, reads: 


where 


(j){z) 


Acoshikoz) + s^(l - smh{kDh)e-^^\^\) iiw>\z\>h 
Acosh^kuz) if ft, > |z| > 0 


A = 
A = 


1 


I S 


sinh{kuw) \ekn 

GZcTla i. U 

eko 


sinh(fcj)ft)e 

efcb 


( 21 ) 


( 22 ) 

(23) 


In order to obtain the modified Smoluchowski equation for the coated capillary one must take into account the fact 
that also the polymer charge contributes to the Laplacian of the electric potential. Thus, we rewrite eq. [T^ as 


'd-^Ux = '^xP + ^PmQux + eEx-^tj) + eZsExUsQ . 


(24) 


The first two terms in the r.h.s. of eq. 24 correspond to those featuring in the r.h.s. of Brinkman equation of porous 
media eq. After substituting the explicit representation eq. [^of (p we recast the Smoluchowski equation as: 


Ux{z) = Acosh(fcpz) -I- ^ if ft > |z| 


Q ^2 

AL 

dz^ 


77 K J < 

■§^Ux{z) - t3^Ux{z) = ^k%Acosh{kDz) 

ez^E^n, ginh(fc^/l)e-fc^l^l 


eZs^^rig+'V j;P 

V 


if IC > |z| > ft 


where we defined the new characteristic inverse length, /3, as 


/3 = 



(25) 
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related to the permeability, k,hy k = 

The solution with the properties that the fluid velocity and its hrst derivative are continuous at z = ±h and vanishes 
at the walls, z = ±rc, can be written in the following form: 

f Ux{z) = Gcosh(fc£)z) + F + if h > |z| > 0 

< U 2 ,(z) = Gcosh(/?z) + Z?cosh(A:£)z) + (26) 

[ itw>\z\>h 

and the constants D, C, G, F, F[ are determined by the set of equations: 

H = -f"^^sinhikDh), 

G= ° 

C cosh(/3'u;) + = —D cosh(fc£iw) — 

< I eZsE^ris+V 

G/3 sinh(/3/i) — Bj3e~^^ = —/cd-D sinh(A:£)/i) + kDe~^^^F[ 

FkuG smA{kjjh) + 

F = G cosh(/3/i) + D coshlkoh) + + Be~^^ 

. -Gcosh(fc,,h) - 


The volumetric flux, neglecting the variations of the densities, is given by Im = Pm dzuz(z) = Pm^ where 

2G G 

<& = -—sinh(fc£)h) + 2F/i + 2—[sinh(/3r(;) — sinh(,5h)] 
ko P 

+ 2^[smh{kD'w) — sinh(fc£ih)] 
kn 

j3 ko 

eZsExTls+VxP , , Vp 3 

- 2- — - {w-h)F—hF 

3p 


Some simple limits can be analyzed: 

• when the polymer charge vanishes and kjyh^ 1 (thin EDL regime) the velocity at midpoint is 

((p ~ cosh(/Z/i) 

pko cosh.(Pw) 

_VpJ_ _ cosh(/3/i) \ 
p \ 2 cosh(/3w) / 

^ if/3h»i 


In the case /3/i <C 1 one recovers the standard EOF plus the pressure induced flow in a slit having the full 
width 2w; in the case I3h^ 1 the EOF is exponentially suppressed at midpoint, while the pressure induced flow 
corresponds to a slit having reduced width 2h. 

If the polymer charge vanishes and kow <C I (thick EDL regime), in the absence of pressure gradient the velocity 
at midpoint turns out to be 

„ (n\ ~ cosh(^\ 

^ pI3'^w\ cosh(/3w)/ 
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which should be contrasted with the polymer-free case results ■Ua;(0) « in the case 

/3/i 3> 1 and /3h ^ 1, respectively. 

The majority of the existing studies [20l |22l [23l |32] have focused on the small Debye length regime (molarities of the 
order IM) , corresponding to very thin EDL, because MD simulations at lower concentrations face severe problems 
of statistical accuracy due to the small number of ions considered. On the other hand, the LBM being based on 
the phase space distribution functions, does not run into such a difficulty and one can access the low concentration 
regime. We shall use the analytical method in the forthcoming section]^ in order to interpret the numerical findings 
and provide a theoretical guide. 


IV. NUMERICAL SET-UP 


The evolution equations for the distribution functions are solved by employing the Lattice Boltzmann method for 
ternary charged mixtures recently proposed by the authors and based on the discretization of the kinetic equations 
on a discrete mesh. A detailed account of the employed numerical methods is provided in ref. |33j so that we refrain 
from repeating the derivation in the present paper for space reasons. The simulated system consists of a slab of fluid 
lattice points (nodes) enclosed by two planes representing the channel walls. The lateral dimensions of the channel, 
expressed in lattice units (l.u.), are 40 x 10 and the channel width, 2w, is 250 l.u.. The presence of grafted polymers is 
represented by two slabs adjacent the two walls each of thickness <5 corresponding to 50 l.u.. Within these two regions 
Ns fixed scatterers, with Ng varying in the range [0,100], occupy at random the 40 x 10 x 50 cells. Flow is in the 
X direction and periodic boundary conditions are used in the x and y directions. The walls are impenetrable to the 
fluid particles and no-slip boundary conditions are implemented by using the bounce back prescription m- 

Throughout this paper we shall use lattice units defined in the following. Lengths are measured in units of lattice 
spacing Ax = 1, the charge e and the mass m are assumed to be unitary, the thermal energy ksT is specified by 
fixing the thermal velocity whose value is vt = s/ksTjm = l/-\/3 in l.u., the kinematic viscosity of water v is set 
equal to 1 l.u., while the dielectric constant is fixed by the Bjerrum length Ib = e^/( 47 refcBlL) and the length 
varies from [10, oo] l.u. as Ng varies between 100 and 0. To simulate the electro-osmotic flow it is necessary to resolve 
the EDL, so that Xu must be sufficiently larger than the lattice spacing. 

In physical units the lattice spacing is Ax^^'^^ = 0.1 nm, Ib = 0.1 nm. = 2.3 nm (corresponding to a 18 mM 

aqueous solution), = 10“®m^s“^ and the mass density of water 18gr/cm^ at room temperature . The physical 

magnitude of the time step is obtained from the relation; 


jyphys 


Atphys 

(Aa;P^!/s)2 


= v 


(31) 


and by using a unit viscosity in lattice units, it follows that = 10 s. 

As in ref. m we tune the parameters of our simulation in order to match those of the experiments where the 
inner walls of a micro-channel composed of a Si 02 glass are coated by Poly(DMA-GMA-MAPS) and a buffer aqueous 
solution of ionic strength of 18 mM made of H 3 PO 4 — NaOH pH 2.5, 6 - e-aminocaproic acid - acetic acid pH 4.4, 
H 3 PO 4 — NaOH pH 7.0, and Bicine-TRIS pH 8.5 fills the pore. The concentration of the charge carriers is obtained 
by imposing the condition of global electro-neutrality for a given surface charge density and with the fixed value of 
the Debye length. 

In order to validate our theoretical estimate of the friction parameter 7 we preliminarily simulated a bulk periodic 
system of size 14 = 10 x 40 x 250 l.u. having the same density, salt concentration and obstacle density as those 
encountered in the channel case. Here the scatterers are distributed uniformly across the whole bulk simulation 
system. We applied a constant uniform force, F, everywhere, in the absence of surface charges and electric fields, and 
measured 7 with the help of the relation between the applied force and the resulting flow velocity according to the 
formula: 7 = The measured 7 compares very well with the theoretical prediction: 


^~^Vg 


(32) 


with kp = 2 l.u., where kp is the strength constant of the drag force introduced in eq. 
l.u. we obtained a value 7 = 10 “^ l.u. (corresponding to 7 P^p® = 10 ^° s“^) and 

pphys ^ QA rim-1). 

Concerning the polymer coating, we choose a polymer slab of thickness d = 5 nm, about twice the Debye length. In 
order to describe charged polymers as experimentally studied by several authors [nnnun] we placed the idealized 
polymers in the regions adjacent to the wall. Following the choice of Monteferrante et al. we varied Ug from 0 


12 By using ris = ^ = 5 x 10 
3 = 10“^ l.u. (corresponding to 


-1-3 





to 5 X 10“^ I.U.. In order to assign such a fractional charge to each scattering center, we desumed the total number 
of polymers coating the walls and the charge of a single polymer from the experiments. For the case of the polymer 
coated channel, the experimentally measured polymer densities ranged between 0.1 and 0.18 gr/cm^. Thus to each of 
the Ns sites randomly fixed in a volume 14 = 40 x 10 x 50 near each wall, we assigned a charge fraction Zg = 10“"^, 
so that the total charge of a polymer layer is Qp = eZsNg. 


V. NUMERICAL RESULTS FOR THE ELECTRO-OSMOTIC CURRENTS 


We first consider the case of neutral polymers: the LBM results for the velocity and potential profiles are reported 
in fig. i and compared to the predictions of the analytical theory of section m In the case of the smaller surface 
charge density S/e = —3.9 x 10“^ l.u., the agreement is better than for S/e = —1.3 x 10“'* l.u. . In fact, the 
linearization of the Poisson-Boltzmann eq. 19 introduces a systematic error which becomes larger when |S| increases, 
as appreciated by comparing the two insets of fig. where the analytical and numerical results for the potential are 
compared. The discrepancy displayed in the right inset has repercussions on the velocity profile, that is appreciably 
lower in the central region in the LBM case. The reason for the discrepancy can be traced back to the small gradient 
assumption used in the constitutive equation for the stress tensor, eq. which ultimately leads to eq. Near 
the walls the velocity gradients are quite large and the Stokes equation might be inadequate to capture the correct 
behavior. The difference between the analytic and the LBM velocity profiles decreases by setting ft, = 65 l.u. in the 
analytical model and ft = 75 l.u. in the LBM. However, such a readjustement is insufficient to produce the same level 
of agreement in the case of larger surface charge, see fig. right panel. Probably one can improve the matching by 
adjusting both the value of E and ft in the analytical model, but we did not pursue further such a program because 
somehow arbitrary, although it could provide a simple and economical tool to scan the overall behavior of the system. 

We consider, now, the case of charged coatings, where the electric potential depends altogether on the surface 
charge, the fixed charges associated with the fixed obstacles and the ionic charges in the EDL. For low values of 
the negative surface charge it is possible to observe a flow reversal in the presence of positive polymer charges. The 
density of mobile charges in the central region |z| < ft has the same sign as the surface charge when E is lower than 
the value 


E* = 


eZgUs 

ko 


(sinh(/c£ift)e — sinh(fc£)U>)e 


_ eZsUs ^kpS 

—eZsUsS 


if koS » 1 
if koS 1, 


(33) 


as one can see from eq. [^and eq. |23| It is worth noting that the overall polymeric charge ZsUgS per unit surface does 
not need to be larger in absolute value than the negative surface charge sitting on the wall in order to induce flow 
reversal in the capillary, consistently with the earlier observation by Hickey et al. [S]. It is important to remark that 
the regime studied in most of the existing literature concerns small values of the Debye length, being typically 

0.2 nm . In this case, the EDL near the wall is not hnely resolved and the comparison with the Smoluchowski theory 
is performed by taking the peak position of the velocity as shear plane, where the velocity assumes its maximum value 
Vsi and its derivative vanishes. In our analytical model conducted at larger values of Ad, instead, we can identify 
such a shear plane with the local maximum of the velocity, that appears in the region ft < |z| < w, as a result of 
the competition between the drag force exerted by the polymers and the electro-osmotic force. Its location is at 
distance A « ~ /3) ln(/3/ftD) from each wall and the velocity decays towards the bulk from such a maximum in 

an exponential fashion. Only for small values of Ad the velocity profile displays two side peaks and no dome at the 
center, fig. displays the results relative to a weakly charged polymer coating, while the remaining parameters are 
identical to those employed in fig. Again, for small surface charge the agreement between the analytical result 
with the rescaled value of ft and the LB simulation is fairly good, but deteriorates at larger values of E for the same 
reasons discussed above. 

As displayed in the inset of fig. left, the potential 4> associated with a small negative surface charge, is non 
monotonic, reflecting the strong inhomogeneity of the ionic charge distribution: the regions adjacent the walls are 
richer of counterions, whereas further away the coions prevail. The direction of the mass flow occurs in the x direction 
in the first region and in the opposite direction in the negatively populated region. Because of the relatively large 
value of the screening length of fig. |^left one can see a large bulge in the center, while for Ad ^ (5 the flow inversion 
takes place only at the walls . 

We obtained the solution of the Stokes eq. [^by imposing the continuity of the derivative of the velocity at z = ±ft. 
However, the latter condition is not necessary in principle, and has been introduced after observing that the LBM 
velocity profiles do not display cusps. Without such a continuity requirement, the analytic velocity profiles obtained 
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by setting the parameter = 0 in eq. display lower values of the velocity in the central region, hence are closer 
to the LB result, but also displays a cusp which is not observed numerically. 

The large kow regime is of particular interest, because in the case of uncoated channels it displays plug-like velocity 
profiles: a very thin region near the walls of thickness Xd where the velocity raises from the zero value at the wall 
to the plateau value Ux{0)- The effect of the obstacles is dramatic: one observes only a single peak near each wall 
and the structure of the velocity profile depends weakly on ku, the relevant length being, now, If the polymer 

coating is sufficiently thick {j36 > 1), besides the peaks near the walls the velocity drops exponentially towards the 
center, with the characteristic length I3~^. The peaks result from the competition between the electroosmotic driving 
force, determining the growth of the velocity from the wall value (respecting the no-slip boundary conditions), and 
the antagonistic frictional force tending to suppress it. 

In this regime the coating appears to be very efficient in suppressing the mass flow, as shown in fig. The 
continuity condition of the derivative of Ux becomes irrelevant in the structure of the analytic solution since the cusp 
is hardly detectable and the difference between the solution with continuous derivative and discontinuous derivative 
is very small. We also remark that in the large ko regime, the flow reversal is characterized by the presence of two 
peaks near the walls, where most of the flow occurs, whereas in the small kjj regime the majority of the flow occurs 
at the center of the slit. 

In fig. [5] we report the normalized mass flow rate as a function of the ratio between the total surface and total 
polymer charge. The theoretical and the simulation results are quite similar, and the difference between the theoretical 
curves h = 75 and h = 65 l.u. is not large. This figure shows that the analytic prediction, as far as channel averaged 
quantity is concerned, is quite accurate and can be used to give quantitative information about the global behavior. 
The right panel of fig. [^displays the mass flow, $, versus the number of charged obstacles, Ng, for different values of 
the surface charge. Clearly, <I> decreases with increasing values of Ng, due to the increased friction. 

In fig. [^we propose a phase diagram using the E, (zgUgS) plane: above the curve the flow is in the positive direction, 
whereas below the curve the flows is reversed. The different curves correspond to different Debye lengths expressed 
in lattice units: Xd = 2.3,10,23,50 (blue, green, red, black). The dashed line represents the curve —S* = ZgUgS the 
electroneutral line, that is, the line obtained by assuming that the charge in the polymer layer equals exactly the 
surface charge. Such a limit is attained in the weak screening regime (Ad > 6 ). For all Debye lengths considered, the 
demarcation curve between direct and reverse flow shows a linear shape. Notice the dependence of the slope of the 
line on the Debye length: the smaller Xd corresponds to steeper lines since the surface charge is screened faster and 
produces a smaller effect. The larger Xd the bigger the polymer charge density where the flow reversal occurs, since 
the effects of the surface extends over larger distances. 

Finally we comment on the charge current, while the mass current is approximately linear with respect to E, the 
electric current is quadratic. Whereas the mass current displays inversion for small values of the surface charge and 
positive adsorbed polymers, the charge current does not. The reason is that the majority of the charge sits near the 
walls where the velocity is also large. The reversal of mass current is due to the fact that this is a flat average of 
Ux{z), while the charge current is an average of Ux(z) weighted by the local charge of carriers, and thus no inversion 
takes place. 


VI. CONCLUSIONS 

In order to investigate the modulation of the electroosmotic flow in polymer coated capillaries we developed a 
phenomenological model and described the grafted polymers by a set of charged fixed obstacles exerting both a 
Coulomb force and a drag force on the fluid species. The electrolytic solution is represented by a ternary mixture 
and its behavior is studied numerically. By analyzing the resulting velocity profiles we found that the coating induces 
features that are not observed in standard EOF, such as non monotonicity, velocity inversion, suppression of the 
plug-like profile in the small Xd regime. This opens the possibility of designing functionalized capillary surfaces in 
order to improve the resolution in capillary electrophoresis. A remarkable feature of the model is that it lends itself 
to an analytical treatment in the case of moderate surface charges, so that the rather complex behavior in terms of a 
reduced set of parameters can be predicted. We have obtained the explicit analytical solution of the linearized model 
and shown that it agrees at semi-quantitative level with the numerical solutions. 

Our approach replaces the complexity of the polymer layer by an assembly of scatterers and thus neglects many 
important aspects such as the deformability of the polymers under the flow, their connectivity, or the possibility of 
forming mushroom or brush structures as discussed by Harden et al. [18] . Experimentally it would be important to 
establish a closer connection between our parameters S and /3 to the degree of polymerization JV (related the polymer 
thickness in the brush regime) and to the number density of grafted polymers onto a flat surface, respectively, instead 
of fixing them by a fitting procedure. Another intriguing aspect that can be investigated by our methods is the 
role played by the charge distribution within the polymer layer on the resulting EOF. As an example. Danger et 


10 


al. [H] used polyelectrolytes of different charge densities to control the EOF, obtained by depositing a first cationic 
polyelectrolyte layer followed by depositing a second polyelectrolyte layer based on anionic copolymer. 

Before concluding, we remark that the differences between the planar and the cylindrical geometry are quantititatve 
rather than qualitative, so that it is possible to extend the present approach to include such a geometry, with an 
analytic treatment slightly more involved and less transparent. 
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FIG. 1: Sketch of the pore system. The near-wall regions represent the polymer coating, the small circles represent the solvent 
particles, the large circles represent counter and co-ions. 
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FIG. 2: Neutral polymer coating. Velocity profiles for /3 = 0.1 l.u. and two values of the surface charge density: E/e = 
— 3.9 X 10“® (upper) and E/e = —1.3 x 10~^ l.u. (lower panel). The simulation curve (circles) is obtained for 2w = 250 and 
h — 75 l.u.. The two theoretical curves are obtained for 2w = 250 l.u. and h = 75 (solid) and 65 l.u. (dashed), respectively. 
Insets: electric potential profiles for each case. 
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FIG. 3: Charged polymer coating. Velocity profiles for p = 0.1 l.u. and E/e = —3.9 x 10“® (npper) and E/e = —1.3 x 10“"^ 
l.u. (lower panel). Symbols as in fig. Insets: electric potential profiles. Notice the flow reversal and the non-monotonic 
character of the potential for small surface charge density (left panel). 
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FIG. 4: Theoretical (black) and simulation (red) velocity profiles, for a slit of width w = 250, polymer thickness 5 = 50 and 
neutral polymer. The curves correspond to /3 = 0 (circles, solid) and j3 = 0.1 l.u. (squares, dashed), and for E/e = —1.3 x 10“"^ 
l.u. and A_d = 7.5 l.u.. 












14 




FIG. 5: Scaled mass flow rate (7o = Ex’ELxLy/v) vs the ratio between the total surface and polymer charge Qp/Qs (upper) 
and vs the number of charged obstacles for different values of the surface charge (lower panel). Data on the left panel correspond 
to simulation (red circles) and theoretical results (black) for h — 75 (continuos line) and 65 l.u. (dashed line), respectively. 
Data on the right panel are from simulations and for Xd = 23 l.u.. Notice the flow reversal for the two lowest values of the 
surface charge. 
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FIG. 6: Phase diagram. The four lines correspond to Ad = 50 (circles), 23 (squares), 10 (diamonds) and 2.3 l.u. (triangles 
up). The EOF is positive above each curve and negative below. The dashed line represents the line —E/e = z^UsS obtained 
when the charge in the polymer layers equals the surface charge (local electroneutrality assumption), a limit valid in the weak 
screening regime (Ad >5). 







